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Abstract 

This is a personal summary of the meeting "Saturation, the Color Glass Condensate 
and Glasma: What Have we Learned from RHIC?" that took place at BNL in May 
2010. The purpose of the meeting was to discuss the status of high density QCD 
and parton saturation, and to review the progress that RHIC has allowed in the 
field. 



1 Introduction 

Understanding nucleus-nucleus collisions at high energies requires a good knowl- 
edge of that part of the nuclear wave functions that describes "low x" degrees 
of freedom. Let us recall indeed that most particles produced in such colli- 
sions have momenta in the GeV range or less, and the partons involved in 
the production process have momentum fraction x = {m±/ ^/s)e^y where ^/s 
is the nucleon-nucleon center of mass energy and y the rapidity of the pro- 
duced particle. Taking m± to be ~1 GeV, one sees that particles produced at 
mid rapidity (y ~ 0) at RHIC probe the nuclear wave functions at values of 
X < 10~^, and at much smaller values of x at forward rapidity. At the full LHC 
energy, values as small as 10^'* will be reached at mid rapidity. These values 
are comparable to the values probed at HERA in deep inelastic collisions on 
protons. 

It is well established that in this regime of small x, the gluons dominate the 
hadron wave functions. The gluon density is large and grows with lowering x, 
but this growth eventually saturates: one reaches then the regime of 'parton 
saturation', which has been much studied over the last decade. Saturation 
leads to a simple structure in the plane (ln(l/a;). In Q^), whose coordinate axis 
are, loosely speaking, the longitudinal and the transverse components of the 
parton momenta: a line separates dense and dilute parton systems. This line 
is parameterized by a 'saturation momentum' Qg, a growing function of 1/x: 
partons with transverse momentum k± ^ Qs are in the dilute regime, those 
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with k± < Qs are in the dense saturated regime. The determination of the 
satm^ation momentum requires the solution of non hnear evolution equations 
that incorporate physics beyond that of the linear evolution equations com- 
monly used to fit deep inelastic scattering data. In the vicinity of saturation, 
perturbation theory breaks down: the large gluon density compensates for the 
weakness of the coupling, making the effective expansion parameter of order 
unity. 

While the gross features of the saturation regime are well identified, the physics 
of saturation, or more generally of dense partonic systems, is a multifaceted 
many-body problem that has been approached using various formalisms, each 
one emphasizing one particular aspect of the problem. These formalisms of- 
ten become rapidly very technical, which makes it difficult to compare their 
physical contents, and relate their predictions to experimental data. In view of 
this, I have chosen in this presentation to put the various facets of parton sat- 
uration, as they emerge from different formalisms, in a balanced perspective. 
I shall keep the discussion simple, my goal being to reveal as much as possible 
the connections between the various points of view in physical terms, rather 
than dwelling on their technical intricacies. I hope the exercise will have some 
usefulness. 

There are excellent reviews on the formal and phenomenological developments 
to which I refer the reader for complete references [1-8] (sec also the relevant 
lectures in [9]). I shall indeed not systematically quote the original literature, 
except in cases where I am using explicitly a specific result. Contributions to 
the workshop will be referred to by the names of their authors. I should how- 
ever mention that, because I have chosen to put emphasis on the conceptual 
issues rather than going through a systematic account of the presentations, I 
shall not be able to cover all the talks that were presented. 

This paper is organized as follows. The next three sections introduce the 
physics of parton saturation from three different perspectives. First I briefly 
review the familiar evolution of parton distributions as a function of x and Q^, 
as describred by linear equations that are solidly rooted in standard techniques 
of perturbative QCD. I also discuss there the early description of parton sat- 
uration as a balance between gluon splitting and gluon recombination. In the 
following section, I introduce further concepts, that grew from the analysis 
of deep inelastic lepton-hadron collisions, namely the eikonal approximation, 
Wilson lines and color dipoles. The study of the propagation of a color dipole 
through the color field of a hadron will provide another point of view on satu- 
ration, which is seen there as arising from the multiple scatterings on the dense 
system of gluons. Finally, I turn to the Color Class Condensate (CCC) where 
the classical fields are given a prominent role. The CCC aims at a complete 
description of the small x part of hadron wave functions that can be used to 
calculate many processes dominated by small x partons. The phenomenologi- 
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cal applications are presented in the next to last section, while the last section 
summarizes the conclusions. 



2 The wave function of a hadron at high energy 



At the beginning of this discussion it is perhaps appropriate to recall that the 
notion of a wave-function for a hadron at high energy suffers from well-known 
ambiguities: it depends on the frame where it is defined, on the gauge chosen, 
with the parton picture emerging more naturally in the light-cone gauge and 
in the infinite momentum frame. A further ambiguity arises often in higher 
order calculations of a given process in the separation of the constituents of 
the hadrons from the probe that is used to measure them. 

The wave function of a hadron is commonly characterized in terms of 'partons' 
carrying momentum k = {kz, kj_), with the longitudinal momentum kz given 
as a fraction x of the momentum of the parent hadron, kz = xPz- We shall 
often in the following refer to light cone coordinates, with = (t ± z)/ ^J2. 
In these coordinates, a right mover parton has momentum /c"*" = xP^ . 




Fig. 1. The elementary branching process: a parton (e.g a valence quark repre- 
sented here by the straight line) emits a gluon (the wavy line) with momentum 
kz = {xPz , k± ) , with Pz the momentum of the parent hadron to which the initial 
parton belongs. 

The basic phenomenon that controls the evolution of the wave functions is the 
branching of partons, with the elementary process displayed in Fig. 1. This 
process, which corresponds to the radiation of a soft gluon from either a quark 
or a gluon, occurs with a probability 

dP^^^^, (1) 



TT 



ki X 



with Cr = Ca = Nc for the radiation from a gluon, and Cr = Cp = {N^ — 
l)/2Nc for a radiation from a quark. As revealed by Eq. (1), this probability 
is enhanced when the emitted gluon carries a small transverse momentum or 
a small energy fraction. 

The state of a hadron is built up from successive sphttings of partons starting 
from the valence quarks. For just a valence quark, leading order perturbation 
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theory yields the integrated gluon distribution xG{x, Q'^) as 



xGM^)^^ln(-^). (2) 

y^QCDj 

Roughly speaking, xG{x, Q^) counts the number of gluons in the hadron wave 
function (here the valence quark) with longitudinal momentum between xP~^ 
and {x + dx)P~^, and localized in the transverse plane to a region of size 
Ax± ~ 1/Q (corresponding to the wavelength of the partons with the highest 
transverse momenta). The infrared cut-off ~ ^qcd accounts for the fact that 
the parton description ceases to make sense for partons that have wavelengths 
larger than the typical confinement scale ro ~ 1/Aqcd- 

It is convenient to relate xG{x, Q"^) to a so-called "unintegrated parton distri- 
bution" (p{x,k±) more directly related to the phase space density. We write 

xG{x,Q') fQ , dN 



with dy — dx/x and 

dN 2{Nl - 1) 



dyd^k^d^b (27r)^ 



ip{x,kA_). (4) 



As is clear from its definition (p{x^ k^) denotes the density in the transverse 
plane of gluons with transverse momentum k± and a definite spin and color. 

The density of partons in the cloud of gluons that surrounds a valence quark 
is not a fixed quantity but a quantity that depends on the resolution with 
which one is probing the gluon cloud, as one can already see from Eq. (2). If 
one increases the resolution, by increasing Q^, on sees more and more partons, 
i.e., the parton density increases. The formula (2) has been obtained by using 
leading order perturbation theory, that is, by taking into account a single 
branching. But even if is small, successive branchings must be taken into 
account when InQ^ becomes large: a large value of is primarily achieved 
through successive emissions of partons with moderate transverse momenta. 
In technical terms, when In becomes of order unity, higher order terms 
become significant and must be resummed. This resummation is achieved by 
the DGLAP equation [10], which we write schematically as 

where P{x/z), the 'splitting function', gives the probability that a daughter 
parton with momentum x is produced by the splitting of a parent parton with 
momentum z. 
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The DGLAP equation leads to an increase of the parton density with in- 
creasing . However, this increase is slow, involving typically InQ^. Since 
the size of the added partons decreases as l/Q^, the area occupied by these 
new partons in the transverse plane eventually decreases with increasing Q'^ . 
Thus, even though the density increases, the system of partons produced by 
the DGLAP evolution is effectively more and more dihitc with the partons 
effectively weakly coupled. As grows (with x kept not too small) perturba- 
tion theory becomes more and more reliable in describing the changes in the 
hadron wave functions. 

When one increases the rapidity (or equivalently the energy), 1/x decreases 
and one eventually reaches a regime where new corrections become important: 
when Q;5ln(l/a;) becomes of order unity, the corresponding large logarithms 
need to be resummed. This new resummation is achieved by the BFKL equa- 
tion [11]. Before we turn to this equation, let us observe that there is a regime 
where both large and small x effects are simultaneously taken into ac- 
count, albeit partially. This is the so-called Double Logarithmic Limit (DLL), 
where ln(l/a;) InQ^ — 1. The evolution equation in this limit (and at fixed 
coupling) reads 

xG'(x,g') = "^xG'(x,g2), (6) 



ain(l/a;)aing2 ^ ' - ^ 7^ 

whose solution is xG{x,Q'^) oc exp ^^2^as ^ Q^-j , with as = cksCa/tt. 

This solution reveals the growth of the structure function at small x, a growth 
which is however milder than that predicted by the full BFKL equation to 
which we now turn. 

Written as an equation for the unintegrated gluon density, the BFKL equation 
takes the form (y — \a(l/x)) 



a. 



dy J 71 pl{k^-p^] 



2 



(7) 



As it should, it contains the DLL obtained when k"^ is sufficiently large. But 
the most remarkable feature of the BFKL evolution is the exponential growth 
that it predicts for the gluon density: 

^{y,kl)r.e^^^^, (8) 

with a; = 4 In 2 (in leading order) . This explosive growth may be traced back 
to important properties of the BFKL cascade: the evolution takes place at 
(approximately) fixed transverse momentum, and the longitudinal momenta 
for the successive gluon emissions i and i + 1 are ordered, namely Xj+i < Xi. 
Since the typical lifetime of a fluctuation is Ati ~ 2xiP~^ /k'j_, Xj+i < Xi implies 
Ati > Aij+i. Thus, the nth gluon is effectively emitted from a strong color 
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source, made from the gluons emitted at Xi^n, which overlap in the transverse 
plane and which can be considered as frozen during this last emission (At„ < 
Atj<„). This enhanced color charge triggers further emission leading to a kind 
of a chain reaction responsible for the exponential increase [12] . 

It was recognized early on that this growth of the gluon density, predicted 
by the linear BFKL equation, could not go on for ever, and various mecha- 
nisms leading to a 'saturation' of the process have been looked for. The early 
approaches to saturation invoked a non linear contribution to the evolution 
equation [13,14] and leads (schematically) to an equation of the form 



92 xG(x,Q2) 
(91n(l/a;)(91nQ2 



as xG{x, Q"^) — 



9 

16' 



_2 ^[xG{x,Q^)f 



(9) 



which differs from Eq. (6) by the second term accounting for 'gluon recombina- 
tion'. (The square of the gluon distribution is an approximation for a 2-point 
gluon distribution, in the form used in [15].) This "kinetic" vision of gluon 
saturation suggests immediately the existence of a characteristic momentum 
scale at which the processes of gluon emission and gluon recombination bal- 
ance each other. This saturation momentum Qs is (parametrically) given by 



2\ xG{x, Qg) 



7tR^ 



(10) 



Referring to Eq. (3), one sees that at saturation the phase space density of 
modes with k\ < is large, of order l/ag. This is an important feature of 
saturation, which we shall often refer to. 



log(^-') 




#- log(e') 



OCD 



Fig. 2. The various regimes in the plane ln{l/x),lnQ'^ (from F. Gelis). At large Q^, 
and moderate values of x, partons are weakly coupled and form a dilute system. 
At small values of x, partons form a dense ('saturated') system and are strongly 
coupled (shaded area). The region with < ^qcd genuinely non perturbative 
and not amenable to a weak coupling description. 



6 



In slightly more general terms, the onset of saturation coincides with a break- 
down of perturbation theory where gluon interaction energies become compa- 
rable to their transverse kinetic energies, that is d'^ ~ as{A^)Q, where {A^)q 
denotes the fluctuations of the gauge fields with transverse momenta up to 
Q, {A^)q ~ xG{x, Q'^^I'k'P?. Thus, as illustrated in Fig. 2, the saturation mo- 
mentum separates partons into dilute modes with ^ Qg that are weakly 
coupled, and modes with < Qs that are strongly coupled because densely 
populated. 



3 Color dipoles, Wilson lines 

Another view of saturation, that will eventually lead us to the most elaborate 
non linear evolution equations, builds on a picture commonly used in the 
analysis of lepton-hadron deep-inelastic scattering (DIS). In an appropriate 
frame, one can describe the interaction of the virtual photon with the hadron 
as the interaction of a color qq dipole (emerging form the photon) with the 
color field of the hadron (see Fig. 3). The factor in the interaction cross section 
that is relevant for the present discussion is o^^^ix^ r^), the total dipole-hadron 
cross-section (that we assume for simplicity here to be a function of x and r_L 
- for a recent and thorough analysis of this picture see [16]). This dipole cross 
section can be calculated in the eikonal approximation, with the size r ^ of the 
dipole remaining unchanged during the interaction. In this approximation, the 
S-matrix for the scattering of a quark moving in the negative z direction is 
given by the Wilson line 



where V denotes an ordering along the x" axis, and is the classical (frozen) 
color field of the hadron moving close to the speed of light in the ^z direction. 
The S-matrix for the scattering of the dipole contains another, complex con- 
jugate, Wilson line. Keeping in mind that after averaging over the field of the 
hadron the S'-matrix will be real, we can write the total dipole cross section 



In fact we shall ignore here the impact parameter dependence and write simply 

S{r It is also customary to set S* = 1 — iV, with denoting the imaginary 
part of the forward scattering amplitude. (The real part is negligible at high 
energy, and would disappear anyway in the gaussian averages that we are 
going to perform.) Obviously, when — > 0, the scattering amplitude vanishes. 




(11) 



as a,,p = 2/^26(1- 5(6, r^)) 



5(6, rx) = ^Tr (c/(6 + ^) lj\h - ^)) 



(12) 
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independently of the field configuration, a characteristic property of the dipole 
interaction often referred to as 'color transparency'. 




Fig. 3. Pictorial view of the interaction of a virtual photon and a hadron. In the cho- 
sen frame, the virtual photon splits into a qq pair, a color dipole, that subsequently 
interacts with the fluctuating color field of the hadron. At high energy, the propaga- 
tion of the dipole through the hadron can be treated in the eikonal approximation, 
i.e., by Wilson lines, and during the passage of the dipole through the hadron, it 
"sees" a frozen color field ^ schematically indicated by the shaded area. 

As a first orientation on the effect of the averaging over the field configura- 
tions (denoted by (■ ■ ■ ) in Eq. (12)), we may assume that their distribution is 
gaussian. Then Wick's theorem leads immediately to the following expression 

-g'Cp j (^^4 — A{k^)\ , (13) 

which involves the correlator of the color charges of the hadron: {p"'{k±)p'^{—k±)) = 
SahT^R^ p\{ki_). Specifically, we have assumed, for the correlator of two fields 
(after integrating over the longitudinal coordinate x", A'^{xj_) = J dx~A'^{x~, x±), 
and assuming that fields at different x~ are not correlated), 

{AU^±)At{yJ) = Sab J ^2^-^^^A{k±), (14) 

with r± = x± — yj^. In the MV model, to be discussed below, all corre- 
lations between the color charges in the transverse plane are ignored. In 
that case p\{k) = /i^ is a constant, p\ = N""""^ with p^ corresponding 
to a single valence quark, for which {pa{k±)pb{—k±)) = dahO^ I'^N^., so that 
/i2 = (7V(27ri?2iV,). 

The momentum integral in Eq. (13) is reminiscent of Coulomb scattering over 
a distribution of charges, with the second term accounting for the interference 
responsible for color transparency. In fact the integral is very well approxi- 
mated by its contribution for < l/r_|_, 

J 27r 2 k\ 167r ^ ' 
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The infrared cut-off vq ~ 1/A.qcd has been introduced to account for the 
screening of the Coulomb field over such distance scales. (We assume rj_ <^ 
throughout.) The factor r]_ in the equation above is characteristic of dipole 
interaction. The total cross section will depend on the ratio between rj_ and 
another length scale, call it r^, that is determined entirely by the dynamics of 
the field with which the dipole interact. More concretely, we may write, 

S{r^) = e-^''-' (16) 

with Qs = I/ts (the factor 4 is conventional). In the case where the field is 
created by a random distribution of (A'"''"' = ANc) valence quarks, we have 

2 

Ql = Ql]n^ = ni^lxGix, 1/rl), Ql = a^Cpf^l, (17) 

where we have absorbed the logarithm into the gluon density of a single valence 
quark at the scale ~ l/r^^ (see Eq. (2)). Equation (16) exhibits the expected 
change of regime in the interaction of the dipole with the field of the hadron. 
A small dipole, with r± <^ Vg, is little affected, and its scattering amplitude 
measures directly the gluon density. A large dipole (r^ ^ r^,) on the other 
hand is strongly absorbed: its cross section saturates to the black disk limit, 
and it is not capable to resolve the parton structure of the field in which it 
propagates. 

It is in fact interesting to push this discussion a bit further, and note that the 
inelastic cross section can be written as ainei = J d'^b (1 — S^{b, r±)), with 

S\r^) = e-^ = e-Mr±)n{b) ^ ^-l/x^ ^^g^ 

where n{h) is the density of valence quarks per unit transverse area at impact 
parameter h (we take it independent of 6, n{h) = N'""''' /irR"^), and the last 
equality expresses as a survival probability with A a mean free path, and 
L the path length. The dipole cross section, corresponding to the scattering 
off a single valence quark, is given by 

ao{r^)^'^rlxG{x,l/rl). (19) 

So what the eikonal approximation does, when coupled to the gaussian av- 
eraging over the field fluctuations, is to express the survival probability as 
the exponential (18), where the cross section ctq is the two gluon exchange 
estimate (19) of the dipole cross section. This allows one to view saturation 
as resulting from the multiple scatterings that the dipole undergoes as it tra- 
verses the hadron. Referring back to the origin of the calculation in terms of 
Wilson lines, one sees that when multiple scatterings are important one can no 
longer expand the Wilson lines: The weak field approximation breaks down. 
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We arc in a regime of strong field, an aspect of saturation that we have aheady 
encountered. 

It may be useful to note that much of the previous discussion would apply as 
well to the scattering of an electric dipole d on a collection of random charges. 
The Hamiltonian of this system would read H — —d ■ E and the associated 
survival probability would be given by 

1_ = ^E^T^) (20) 

where T denotes the time the dipole is exposed to the field (T plays the role 
of the length of the x" integration). In each passage of the dipole through 
the field, the field is constant so that the wave function of the dipole is just 
modified by a phase (1 assume that d and E are parallel for simplicity). But 
if we take the field to be a random variable, and take the distribution of this 
variable to be a gaussian, we arrive at Eq. (20). To make the connection with 
the previous calculation involving color fields, we note that the fluctuation of 
the color electric field = = VjA"*" can be obtained from the gluon 
distribution (again for the field created by a random distribution of valence 
quarks) as (using Eq. (14)) 

{Eix^)-E{x^)) = ^^xG{x,Q'), (21) 
so that (7o can indeed be written as 

r^ib)<^o^~{E{x^)-E{x^))rl. (22) 



This picture of the dipole as a probe of the hadron wave function leads one 
to expect a rather simple property that will be more sharply defined later. 
We note that all the dependence on the wave function is contained in the 
scale Tg that characterizes the change of regime between the dilute regime 
and the saturated one. The energy dependence of is determined by the 
dynamics of gauge fields, but we expect to appear in the cross section 
only in the ratio r^/r^. In fact, a very simple model for the dipole cross- 
section has been proposed by Golec-Biernat and Wiisthoff (GBW) [17] with 
the parameterization 

<7dip(^,r-x) = a[l-e-i«^(-H]. (23) 

Quite remarkably, with Q^x) = Qq (xq/x)^, this leads to a very good descrip- 
tion of HERA data at x < 10^^ and moderate [18]. The fact that the x 
dependence of the cross section enters only in the definition of the basic scale 
Ts has been dubbed 'geometrical scaling'. 
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Fig. 4. Basic effect of the evolution on the dipole. An increase of rapidity opens the 
phase space for the emission of one gluon. This contributes to the screening of the 
charge, modifying the effective interaction of the dipole as it traverses the color field 
of the nucleus. In the large Nc limit, the dipole-gluon system is equivalent to a pair 
of elementary dipoles, with endpoints a;, z and z,y. 

Leaving aside phenomenology, we note that r^, or equivalently the saturation 
momentum Qs, is determined by the dynamics of the gauge fields, i.e., from 
QCD. In particular, the energy dependence of Qs is expected to follow from 
the non linear evolution equations that capture the phenomenon of saturation. 
These equations exist in several versions. The dipole picture provides perhaps 
the most direct and intuitive approach to these equations. In this picture, the 
main agent responsible for the evolution is the emission of a gluon by a color 
dipole, which takes place when the dipole is boosted to a high enough rapidity. 
This is depicted in Fig. 4. After the emission of a gluon, the original dipole 
turns into a dipole-gluon system, whose propagation in the field of the hadron 
differs from that of the original dipole. The evolution equation is precisely the 
equation that accounts for this modified propagation. For the correlator of two 
Wilson lines, this equation takes the form 

9y(Tr {ulUy)), = -^2/ d'^l^^y^ (NcTr (f/tf/,) - Tr (f/tf/.) Tr (f/tf/,)), 

(24) 

where JC^yz dPz dY, with 

(33 — y)'^ 

= {x-zy{y-zy' ^^^^ 

is the coordinate space version of the branching probability (1). This equation 
is in fact the first in an infinite hierarchy of equations for the correlators of 
an arbitrary number of Wilson lines that was obtained by Balitsky [19]. The 
JIMWLK equation [21] is a functional equation that describes the evolution of 
the generating functional of all these correlators. Its physical content is iden- 
tical to that of Balitsky's hierarchy. The Balitsky-Kovchegov (BK) equation 
[20] is obtained by assuming a factorization of the non-linear term 

(Tr (ulUz) Tr {ulUy)) ^ (Tr (f/It/.))(Tr (f/tf/,)), (26) 

which is usually justified on the basis of large N^. arguments. This factorization 
allows to close the Balitsky hierarchy at the level of the 2-point function. The 
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most elaborate calculations to date involve solving the BK equation, with 
running coupling (see the contributions by Kovchegov and by Albacete). For 
numerical solution of JIMWLK and comparison with BK see the discussion 
in [2] and references therein. 

We have presented the evolution equations as reflecting the change in the 
internal structure of the dipole that probe the hadron wavefunction under a 
boost. Of course, moving to a frame where the dipole stays intact and the 
boost affects only the hadron, we can view the emission of the extra gluon 
as part of the modification of the hadron wave function. Thus the evolution 
equations can be as well interpreted as reflecting the evolution of the hadron 
wavefunctions with increasing energy. 




Fig. 5. The dipole scattering amplitude N{r) as a function of the dipole size for 
various rapidities y(from Ref. [25]). The curves are obtained by solving the BK 
equation with fixed coupling (full lines) or running coupling (dashed lines). The 
travelling waves are clearly visible, with a 'saturation front' joining small and large 
dipoles at a scale rs{Y) that decreases with Y. Note that the evolution is slower for 
running coupling than for fixed coupling. 

With the change of variable = 1 — S*, we can write the BK equation as an 
equation for the scattering amplitude A^: 

dyN^y = — / ^ JCxyz {N^z + N,y - N^y - N^yN.y) . (27) 

TT J ZTT 

Note that the first three terms (those linear in A^) yield the BFKL equation, 
in a form that underlines the structure of the hadron wave function in terms 
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of color dipoles [22]. Interesting connections have been established between 
this equation and equations that occur in problems of classical statistical me- 
chanics, e.g. in the study of reaction-diffusion processes (see the review by 
Munier [8]). In particular, it was shown that the asymptotic solutions of the 
BK equation have an universal behavior, in the form of traveling waves [23] 
that represent the evolution of the saturation front with rapidity (see Fig. 5 
for an illustration) : 

N{r,Y) = N{r -rs{Y)). (28) 
The property of geometric scaling naturally emerges from such solutions. 

4 Classical fields, CGC 

We have seen classical fields emerging at several occasions in the preceding 
discussion. For instance, we often used implicitly an adiabatic approximation 
in which during the time of interaction between a dipole and a hadron, the 
color field of the hadron can be considered as frozen, i.e., its fluctuations 
can be ignored. Another important feature has emerged: near saturation the 
color fields are 'strong', with amplitude A ~ 1/g, and they provide a natural 
description of highly occupied, strongly coupled, modes. 

The Color Glass Condensate formalism puts these classical fields at the heart 
of all considerations. Its goal is to provide a complete description of the small 
X component of the nuclear wave functions, thereby allowing, among other 
things, the calculation of observables that control the early stages of nucleus- 
nucleus collisions. 

A key idea in the approach is the separation of the small x degrees of freedom, 
treated as classical fields, from the sources that give rise to these fields and 
that are frozen during the collision. One may imagine the various field con- 
figurations at a given rapidity Y to be described by a wave-function $y[A]. 
Because the field does not change during the collision, the average over the 
field configurations will naturally involve the square of the wave function, i.e., 
a probability (at least for sufficiently inclusive quantities). Schematically, the 
calculation of observables at some rapidity Y will then take the form 

{■■■)y = JvA\<i>y[A]mA\...\A). (29) 

But how can we determine |$y[74]|^? 

A daring step towards answering this question was taken by McLerran and 
Venugopalan. Their model implements the separation of scale in the most 
economical way, by taking the field A^ as the solution of the classical Yang- 
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Mills equations: 

[D^,F>^'']^r, (30) 

in the presence of a frozen color source . Typically, J'^{x~, x±) = S^~^p{x~, x±) 
for a nucleus moving at nearly the speed of light in the positive direction. 
The question of determining the distribution of field configurations |$y[y4]p 
is reduced to that of finding the distribution of color charges p, traditionally 
denoted by W[p\. In the MV model, this is taken to be a gaussian distribution, 
characterized by the density-density correlation function 

{p'^{x-,x^)p\y-,y^)) = 6abS{x~-y-)S^'\x^-y^)p\x-). (31) 

That is, one assumes that color sources are completely un correlated. Argu- 
ments that motivate such an ansatz can be formulated in the case of a heavy 
nucleus [24]. 

The gluon distribution provided by the MV model can be easily calculated. It 
provides a non-abelian generalization of the Weizsacker- Williams (WW) field 
of electrically charged particles. This gluon distribution is closely related to 
the correlator S{r±) of two Wilson lines, in the adjoint representation. This 
correlator is given by Eq. (13) in which one replaces Cp by Ca — Nc, with, 
accordingly 

Q's^j;^n{b)xG{x,Ql). (32) 
The actual WW gluon distribution, ipA, is given by 

(33) 

with N{k±) the Fourier transform of N{r±) = 1 — S{r±). Note that this 
distribution is independent of x (since each valence quark produces an x- 
independent density). The behaviors for small and large transverse momenta 
are interesting. We have 

1 n2 1 Q2 

Mk± < Qs) ~ —iTf In ^, Mk± > Qs) ~ -^ji- (34) 

At small momentum, (pA{ki. ^ Qs) exhibits saturation. It is in fact believed 
that the small momentum behavior in Eq. (34) is generic, and indeed, as we 
shall see shortly it is essentially preserved during the evolution. It is quite 
remarkable that this feature emerges from the non abehan WW field calcu- 
lated with uncorrelated charges. Note that in this small momentum regime, 
if A ~ l/^s, reflecting the large occupation of the low momentum modes, as 
expected at saturation. At large momenta, the distribution decays as Ql/k^, 
in agreement with perturbation theory and the additivity of the gluon distri- 
butions of the various sources. 
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Fig. 6. The unintegrated (WW) gluon distribution obtained as solution of the 
BK equation (from Ref. [26]). Initial conditions are GBW (solid lines), MV with 

q2 = 4GeV2 (dashed lines) and MV with Ql = 100 GeV^ (dotted lines). The three 
sets of curves correspond the rapidities y = 0, 5, 10, respectively, from left to right. 



The gluon distribution obtained from the MV model is often used as initial 
condition in the equations describing the evolution towards small x. As an il- 
lustration, we show in Fig. 6 solutions of the BK equation that use MV gluon 
distribution as initial input. The distribution is nearly constant (in a logarith- 
mic scale) up to the knee that corresponds to the saturation momentum which 
grows with increasing rapidity. The evolution produces a mild increases of the 
height of the plateau at small k_[_ , but the dominant effect is the increase in the 
width of this plateau that reflects the increase in the saturation momentum. 
An important feature that is not made obvious by the plot concerns the mod- 
ification of the power law at large momenta, which is of the form {Qs/k±)'^'^ , 
with 7 ?s 0.63. The numerical calculations reveal that this modification of the 
spectrum extends to fairly large values of the transverse momentum. This is 
an intriguing feature, which seems to imply that the approach to the standard 
perturbative QCD behavior is delayed to unexpectedly large momenta. 

The WW gluon distribution is not the distribution that enters naturally the 
calculation of many observables. For instance the spectrum of gluons produced 
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in a nucleus-nucleus collisions is given by 

This factorized expression (strictly valid only in the case where one of the 
two colliding systems is dilute) expresses the result in terms of an elementary 
cross section ~ (y.s/p\ and an integral over unintegrated distribution functions 
defined as 

M^.k^) = 1^ / rf'a;xe-^^-^-(Trt/t(0)t/(a^x))y. (36) 
There is a simple relation between 0a (a^, k±) and (Pa{x, k±) (see Eqs. (33)): 

M^, k±) kl^l^Aix, k±). (37) 

Note that the total multiplicity obtained by integrating the spectrum (35) over 
Pj^ is plagued with a logarithmic infrared divergence that is usually regulated 
by giving a small mass to the emitted gluon (or introducing a cut-off). 

The assumption of uncorrelated charges of the MV model is a crude approxi- 
mation: the evolution rapidly builds up correlations between the color charges. 
It is interesting to analyze how the various stages of the evolution modify the 
density-density correlation function ^^{k) (sec Eq. (14)). In the MV model, 
at the start of the evolution, /i^(fe) = /^^ constant. In the BFKL regime, 
near saturation, /x^(fc) ~ \J ^\k'^e^°''^ . Finally in the deep saturation region, 
li{k) ~ k^{Y — Ys{k)). A physical discussion of this evolution of the charge 
correlation function in terms of the gradual build up of the screening of color 
charges is presented in Ref. [27]. 

At this point we shall introduce the second important ingredient in the CGC 
formalism. It is based on the observation that the separation between charges 
and fields, that is exploited in the MV model, involves a dividing scale between 
degrees of freedom: the field describes degrees of freedom with some particular 
value of X, while the color charge, and their correlations, are determined from 
degrees of freedom with x- values > x. As we move to lower and lower values 
of X more and more degrees of freedom are treated as random charges (cf. the 
discussion of the BFKL cascade earlier), and the corresponding correlators are 
modified. The requirement that observables should remain independent of the 
dividing scale can be implemented (to leading logarithm accuracy) as a renor- 
malization group equation for the distribution PFy[p]. This renormalization 
group equation is nothing but the JIMWLK equation. 

A further progress concerns the factorization theorem presented by Gelis. To 
keep the discussion concrete, let me focus on the calculation of the energy 
momentum tensor at the early stage of a nucleus-nucleus collision. In the 
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CGC formalism, this is given by (r is the proper time and 77 the space-time 
rapidity) 



Pi 



Wo 



P2 



(38) 



where T^^ denote the contribution calculated from the classical fields gener- 
ated by the distributions of color charges pi and p2 of the respective nuclei. 
Equation (38) exhibits the factorization alluded to: the distributions W[p] 
play a role analogous to the parton distributions in the collinear factorization 
scheme; they determine the classical field which enters the calculation of T^^. 
Because the weight functions W[p\ obey the JIMWLK equations, this formula 
resums all contributions of the form (a^. ln(l/x))". Besides, the dynamics of 
the classical field is treated exactly, which is necessary when the sources are 
strong, p ~ 1/g. Note that, in this approach, correlations, and fluctuations, 
are introduced through the averaging over the color charges. 



5 Phenomenology 



We arc now ready to examine some of the phenomena observed at RHIC, 
and discuss to which extent these may reveal features of saturation discussed 
in the previous sections. Two important aspects need to be kept in mind in 
going through this discussion. The first one is that the phenomenology is in 
most cases blind to many of the details of the theory, and tests will therefore 
concern mainly general qualitative trends. The second aspect is that the solu- 
tion of the most sophisticated evolution equation is technically difficult, and 
approximations are often used that tend to obscure the interpretation. 



5.1 The basic ingredients 



Much of the phenomenology based on the saturation picture is driven by 
the saturation scale and its dependence on energy and impact parameter. 
TypicaUy, 

Ql-Ql(-)\ Ql{h) ^ Ql{Q)Tj,{bl (39) 

where 7^(6) = / dzn{b,z), with n{b,z) the nucleon density. These are of 
course only approximate relations, and the energy dependence of Qs, for in- 
stance, can be obtained by solving an evolution equation such as the BK 
equation, or from the more sophisticated treatment presented in Ref. [28]. 
The treatment of the impact parameter remains in most cases rather crude, 
and even the simple additive geometrical approximation suggested by the sec- 
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ond formula above is rarely implemented. In many cases one assumes that 
nuclei are flat uniform cylinders, for which ^^(6) = A/ttR"^. 

A further important ingredient is the running of the couphng constant. In the 
simplest treatments, one usually chooses the scale of the running coupling to be 
the saturation momentum itself, i.e., as — as{Qs)- But there are now sophis- 
ticated implementations of running coupling effects in the evolution equations 
(see the contributions by Kovchegov and by Albacete). A generic effect of the 
running of the coupling is to slow down the effects of the evolution. 

The unintegrated gluon distribution functions that enter the calculation of 
most observables have the typical shapes of the functions displayed in Fig. 6. 
A crude approximation that is often used consists in assuming the simple form 
(p{k±) ~ l/tts for k± <^ Qs{x) and (p{k±) ~ Ql{x)/k'j_ for k± Qs{x). When 
multiplied by the phase space factor k\ (see e.g. Eq. (37)), these unintegrated 
distribution functions are peaked at transverse momenta of the order of Qs, 
so that many processes are dominated by transverse momenta of the order of 

Qs. 

Most calculations rely on the so-called /cr-factorization. While this approxi- 
mation can be justified for the spectrum of particles produced in the collision 
of a dilute system on a dense system (it exploits a linearization with respect 
to the color source of the dilute system), it does not hold in the case of the 
collision of two dense, saturated objects, Uke a nucleus-nucleus collision. Re- 
cent numerical studies within the MV model show that the violation of kr 
factorization can be large in the region of small k±_ [36]. In practice, a cut-off 
is used to control this low k±_ region, and this introduces uncertainties which 
arc difficult to quantify (although one may argue that the effect of the cut-off 
may vary little with energy). 

When discussing some particular situations, like for instance the physics in the 
fragmentation regions, one often mixes factorization schemes, in order to deal 
with systems with vastly different gluon densities: one of the colliding nuclei 
is dense and describable by the CGC, the other is dilute, dominated by large 
X partons, and for it coUinear factorization, with standard parton distribution 
functions, seems appropriate. This mixing of schemes, while 'reasonable', has 
not received a deep theoretical justification. 

An important feature of high density gluon systems is contained in various 
correlations. In the CGC approach, most of these correlations are generated 
by the average over the charge distributions. Depending on the process exam- 
ined, various n-point functions need to be calculated. For instance the pro- 
duction of quark-antiquark pairs in pA collisions involves in addition to the 
2 point function, correlators of 3 and 4 Wilson lines, in the fundamental or 
adjoint representations [29] . The same thing happens in the study of di-hadron 
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production in DIS or proton-nucleus collisions [30,31]. The evolution of these 
correlators with rapidity is usually obtained through the BK equation, that is, 
only 2-point correlators are evolved. This means that n-point correlators have 
to be factorized first. There are numerical indications that the full JIMWLK 
equation and the BK equation lead to close results (see e.g. the discussion in 
[2]), but this may not hold for all situations. 

5.2 Deep inelastic scattering and geometrical scaling 

Before we start reviewing RHIC results, let us recall that most of what we know 
about parton densities comes from lepton-hadon DIS as measured in particular 
at HERA. The fits that determine these parton densities are traditionally 
performed within the DGLAP evolution scheme. The state of the art was 
reviewed by Forte. One important issue in the present context is the role and 
importance of small x effects, beyond those included in the initial conditions. 
Calculations at NLO and at NNLO seem capable to account for the data. 
However Forte presented some evidence for tensions in fits that may suggest a 
stronger small x evolution than what is presently taken into account. Another 
issue is whether the data provide any evidence for BFKL evolution. It was 
argued by Kowalski that one can get a good fit to the data using BFKL, once 
NLO corrections are included, in particular through the dependence of the 
power law that controls the growth of the gluon density. 

A remarkable property of the data, which is beautifully captured by the GBW 
model, is the geometrical scaling, with the x-dependence entirely hidden in 
the energy dependence of the saturation scale Qs, as we discussed carher. The 
physical origin of the scaling is still under debate. In the dipole picture, it 
emerges very naturally (to within logarithmic corrections), as we have seen. 
It can also be understood as a universal property of the small x evolution 
equations in their asymptotic (large rapidity) regime [32]. That this property 
emerges also from DGLAP evolution is a priori less obvious, but it is so, as 
argued by Forte: not only DGLAP preserves geometrical scaling if this is built 
in the initial condition, but it is capable of generating scaling through evolution 
for not too small Q^, ^ 5 GeV^. (For a recent discussion see [33].) Thus, 
even though geometrical scaling is very suggestive of the underlying physics 
of saturation, the situation remains somewhat inconclusive. 



5.3 Particle multiplicities in hadronic collisions 

By assuming that the bulk of produced particles comes from the saturated part 
of the wave functions, one immediately gets the following generic behavior for 
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the multiplicity density [34] 

7ri?2 dy " a,(Q2)' ^ > 

with all dependences on energy or centrality being entirely determined by the 
corresponding dependences of Qg- 

Detailed calculations can of course be made beyond this crude estimate. The 
most sophisticated calculations were presented by Albacete. His calculation 
involves solving the BK equation, with running coupling corrections. The en- 
ergy dependence of the saturation scale Qs{x) is determined from the BK 
equation, with parameters (entering the initial conditions) obtained by fitting 
RHIC data (and are compatible also with HERA data). A mixed scheme is 
used to treat the tails of the multiplicity distributions, dominated by large x 
partons. Because of the poor treatment of the impact parameter, there is a 
large uncertainty in the overall normalization, which however depends only 
mildly on energy, so that the prediction for the LHC should be reasonably 
accurate. Another source of uncertainty is the kx factorization on which the 
calculation is based. 

Note that the CGC only predicts the distribution of 'initial gluons', set free 
typically at a proper time r ~ Qj^- Between this early stage and the freeze- 
out, the system undergoes several non-trivial steps: kinetic and chemical equi- 
libration (possibly with additional parton production), hadronization, etc. The 
fact that the CGC multiplicity accounts well for the data seems to imply that 
there is little room left for the late stages of the coUision to contribute signif- 
icantly to the multiplicity. 



5.4 Limiting fragmentation 



By 'limiting fragmentation', one refers to the property of the rapidity distri- 
bution dN/dy to be independent of the collision energy at high energy, in the 
vicinity of the beam rapidity Y — ln(y^/m), i.e., in the fragmentation region. 
The interval in rapidity over which this phenomenon is observed grows with 
energy. Evidence for such a behavior has indeed been found at RHIC and was 
reviewed by W. Busza. 

In the CGC picture one assumes that the phenomenon finds its origin in the 
initial gluon production, and is little affected by subsequent evolution of the 

system. The nucleus that sits in one fragmentation region is a dilute partonic 
system while the other nucleus is in a saturated state. The rise in multiplicity 
when one moves away from the beam rapidity is then attributed to the growth 
of the parton distribution in the dilute nucleus. More quantitatively, this may 
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be seen by using the A;^ factorized formula (35), recalling that the unintegrated 
gluon distribution functions are peaked around their respective Qs- The rele- 
vant partons in the projectile (with y ~ Ybeam) have large Xi = (px/m)e^~^'"='''", 
and hence small Qs, while those in the target have small X2 = (pr/"^)e~^~^'"^'"" 
and hence a large Qs- It follows that the typical transverse momenta of par- 
tons in the projectile are much smaller than those of partons in the target. 
This remark can be exploited to perform the integration over in Eq. (35)) 
and obtain the multiplicity in the form 

^ocxiG(xi,Q^). (41) 
dy 

As anticipated, this depends only on y — Yieam, and weakly on the quantities 
that are involved in fixing the scale Q^. A detailed discussion of this phe- 
nomenon was presented by A. Stasto. 



5. 5 Initial conditions for hydrodynamical evolution 



Simple properties of the initial, saturated, wave functions have interesting con- 
sequences for the initial conditions to be used in hydrodynamical calculations, 
as reviewed by Dumitru. The discussion is again based on the estimate of 
the multiplicity density using the factorized formula (35). Since the mo- 
mentum integration is dominated by transverse momenta smaller than the 
smallest of the saturation momenta of the target and projectile, the result of 
this integration yields the local density of produced gluons in the form 

^«min{g^„g,%}. (42) 

Now each Qs depends (approximately additively) on the number of partici- 
pants in the corresponding nucleus, that is 

QI,<xTa{s), Ql^^Tsis-b). (43) 

This result may be contrasted to that obtained in a standard Glauber calcula- 
tion, where the local density would be simply Ta{s) + Tb{s — b). It follows in 
particular that the eccentricity predicted by the CGC initial conditions is typ- 
ically larger than that obtained with Glauber initial conditions, which affects 
in particular elliptic flow calculations. 



5.6 Long range rapidity correlations 



Long range rapidity correlations can only be produced at very short time. 
This follows from a simple causahty argument recalled by Gehs. The CGC 
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provides a natural explanation for the origin of these correlations (and an 
estimate for the rapidity range over which they should be expected): the color 
field produced at early times (the 'glasma') is arranged in coherent flux tubes 
of radius of order Q~^, and is correlated over a rapidity interval of order aj^. 
However, these correlations are not uniquely characteristic of the CGC, since 
most other models of particle production encompass such long range rapidity 
correlations. This is so, in particular, of the string model, as recalled in Pajares' 
contribution. 

Long range rapidity correlations manifest themselves in a spectacular fash- 
ion in phenomena involving also strong azimuthal correlations. This is the 
so-called 'ridge' phenomenon. In order to get the collimation in azimuthal an- 
gle, collective flow seems to be required, as discussed by Gavin, Shuryak and 
Grassi. A noteworthy contribution was presented by Alver who introduced the 
concept of 'triangular flow', a simple flow effect resulting from fluctuations in 
the initial density distribution, and which had been apparently overlooked in 
all previous analyses. 

Returning to correlations, let us recall that in the CGC framework, most of 
the correlations arise form the averaging over the color sources. As discussed 
by Lappi, when these sources are strong, of order 1/(7, a new ordering in var- 
ious contributions emerges, which differs from the natural ordering in powers 
of the coupling constant. Interestingly, the dominant contribution leads to the 
negative binomial distribution, a well-known feature of multiplicity distribu- 
tions of high energy hadronic coUisions. The parameter k which characterizes 
this negative binomial distribution is found to be a growing function of Qg- 



5.7 Forward rapidity 

At forward rapidity, as we have already discussed, the wave function of the 
projectile is dominated by large x partons, while the target is a dense, sat- 
urated system. One of the early indications of a rapid evolution of the wave 
functions (albeit mostly in the dilute system) were provided by BRAHMS data 
on d-AU coUisions, and the rapid disappearance of the Cronin peak with in- 
creasing rapidity. Both the existence of the Cronin peak and its disappearance 
with increasing rapidity are natural consequences of the CGC. As shown by 
Albacete, BRAHMS spectra can now be understood quantitatively from the 
CGC framework, using BK equation, with parameters adjusted on DIS data. 
However since other competing explanations cannot be completely ruled out 
(see for instance the talk by Strikman), the issue of whether the BRAHMS 
effect constitutes a signature of saturation remains disputable. 

Perhaps more conclusive evidence will come from the study of di-hadron pro- 
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Fig. 7. The absence of the 'away side' jet observed by the STAR experiment (taken 
from the contribution by Braidot). The (blue) points represent the results of the 
calculation by Marquet and Albacete. 

duction in d-Au collisions, as discussed by Marquet. This is indeed a very 
interesting situation where one can probe the very small x component of the 
nucleus {xAy/s = fcie~^^ + ^26"^^ with yi and y2 the rapidities of the produced 
hadrons). The forward double inclusive pion production has been calculated 
using a mixed formalism in which the wave function of the projectile is de- 
scribed by a standard parton distribution function while that of the nucleus 
is described by the CGC (using BK evolution with running coupling). The 
physics one expects is somewhat similar to that of the Cronin effect, namely 
multiple scattering in the dense gluon system. The multiple scattering of the 
leading quarks going through the nucleus is expected to wash out the back- 
to-back correlations between the produced hadrons, leading eventually to the 
disappearance of the away side jet. Such an effect has indeed been observed 
by STAR, and has been reported by Braidot at this meeting (see Fig. 7). 

There are uncertainties in the calculation that need to be clarified before a 
definite conclusion can be drawn. In particular, the calculation involves specific 
correlators of Wilson lines that are at the moment approximated in terms of 
2-point functions (the ones that one can calculate using the BK equation). It 
is unknown whether the neglected higher point correlations are important or 
not. But the phenomenon is very suggestive of an initial state effect that can 
finds its natural explanation in terms of the large gluon density of the nucleus. 
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6 Summciry 



Saturation is a generic property of QCD in the regime of high parton densities. 
The detailed microscopic mechanisms at work in this regime still need to be 
further analyzed to reach a complete understanding, but the gross features are 
well identified. Non linear evolution equations can be derived from QCD by 
following different routes, all based essentially on weak coupling approaches, 
the non pertubative aspects of saturation arising from the large density of 
partons. (Saturation may not be limited to weak coupling, however: as shown 
by lancu, an analogous phenomenon also occurs in strongly coupled gauge 
plasmas, as described by the AdS/CFT correspondence.) The progress on this 
subject over the last decade has been truly impressive, both theoretically and 
experimentally. 

The progress in theory has been in fact continuous, with major steps forward 
at the beginning of the decade, and before, when the non linear evolution 
equations were established. A better control of the saturation scale has been 
obtained, with more accurate determinations of its dependence on the energy 
and the system size. Next to leading order corrections (in particular running 
coupling constant corrections) have been worked out, and when these have 
been implemented carefully, it has led to a better description of the data. 
In this general landscape, the CGC (in the loose sense that it has acquired 
over the years), has established itself as a reference. It has become a use- 
ful organizing principle, suggesting new ways to look at the data and new 
measurements. It also appears as an essential step in building a space-time 
picture of nucleus-nucleus collisions, and this has triggered new theoretical 
developments. In particular, new factorization theorems have been obtained, 
that allow for better controlled calculations of the initial stages of heavy ion 
collisions. The CGC relies on a separation of degrees of freedom into color 
charge and color fields, and on a renormalization group equation that copes 
with the arbitrariness in the choice of the scale at which that separation is 
implemented. This confers the CGC features of an effective theory. One should 
however recognize that the microscopic definition of color charges, and their 
correlations, still hide subtleties that need to be better understood. 

The CGC, or related approaches, have led to a systematic and successful phe- 
nomenology, based on a few basic ingredients: the saturation momentum and 
its variation with energy, size of the systems, and simple (but not always ac- 
curate) approximations, such as the /cT-factorization. Some features of this 
phenomenology are common to other models (and hence not discriminant), 
but one may argue in most cases that the CGC provides a better connection 
to QCD, and the overall picture it provides is more systematic. Thus, although 
one may, rightly, argue that the data do not provide sharp and direct evidence 
for neither the saturation regime, nor the approach to saturation, it is fair to 
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say that HERA and RHIC data lead to a very coherent picture, with several 
phenomena finding their natural interpretation in terms of high density glu- 
onic systems. This is the case in particular of the physics at forward rapidity 
where the smallest x components of the wave functions are being probed at 
RHIC energy. In fact, the beautiful results on di-hadron production that were 
presented at the meeting may well constitute the first, rather direct, evidence 
for a large gluon density effect in nucleus-nucleus coUisons. 

Acknowledgements: I am grateful to J. Albacete, E. lancu, F. Gelis, A. 
Mueller and J.-Y. OUitrault for many useful discussions, and comments on 
the manuscript. I also wish to thank L. McLerran for his patience in waiting 
for this write up. This manuscript was partly completed during my stay as 
Hans. D. Jensen Professor at the Insitute for Theoretical Physics at Heidelberg 
University. 

References 



[I] E. lancu, R. Venugopalan, In *Hwa, R.C. (ed.) et al.: Quark gluon plasma* 
249-3363 . [hep-ph/0303204] . 

[2] H. Weigert, Prog. Part. Nucl. Phys. 55, 461 (2005) [arXiv:hep-ph/0501087]. 

[3] D. N. Triantafyllopoulos, Acta Phys. Polon. B 36 (2005) 3593 [arXiv:hep- 
ph/0511226]. 

[4] J. JaliUan-Marian and Y. V. Kovchegov, Prog. Part. Nucl. Phys. 56, 104 (2006) 
[arXiv:hep-ph/0505052]. 

[5] F. Gelis, E. lancu, J. Jalilian-Marian and R. Venugopalan, arXiv:1002.0333 
[hep-ph]. 

[6] T. Lappi, arXiv: 1003. 1852 [hep-ph]. 

[7] L. Frankfurt, M. Strikman and C. Weiss, Ann. Rev. Nucl. Part. Sci. 55 (2005) 
403 [arXiv:hep-ph/0507286]. 

[8] S. Munier, Phys. Rept. 473 (2009) 1 [arXiv:0901.2823 [hep-ph]]. 

[9] "QCD perspectives on hot and dense matter" . Proceedings, NATO Advanced 
Study Institute, Summer School, Cargese, France, August 6-18, 2001," Kluwer 
Academic Publishers, 2002, Eds. J. P. Blaizot and E. lancu. 

[10] G. AltareUi, G. Parisi, Nucl. Phys. B 126, 298 (1977); V.N. Gribov, L.N. 
Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); ihid. 675 (1972); Yu. Dokshitzer, 
Sov. Phys. JETP 46, 641 (1977). 

[II] L.N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976); E.A. Kuraev, L.N. Lipatov, 
V.S. Fadin, Sov. Phys. JETP 45, 199 (1977); I. Balitsky, L.N. Lipatov, Sov. 
J. Nucl. Phys. 28, 822 (1978). 



25 



[12] A. H. Mueller, "Deep inelastic scattering and small X physics," Lectures given 
at NATO Advanced Study Institute: Frontiers in Particle Physics, Cargese, 
France, 1-13 Aug 1994. 

[13] L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Kept. 100 (1983) 1. 

[14] A. H. Mueller and J. w. Qiu, Nucl. Phys. B 268, 427 (1986). 

[15] J. P. Blaizot and A. H. Mueller, Nucl. Phys. B 289, 847 (1987). 

[16] C. Ewerz, O. Nachtmann, Annals Phys. 322 (2007) 1635-1669, [hep- 
ph/0404254]; ibid. 1670-1726, [hep-ph/0604087]. 

[17] K. J. Golec-Biernat and M. Wusthoff, Phys. Rev. D 59, 014017 (1998) 
[arXiv:hep-ph/9807513]; Phys. Rev. D 60, 114023 (1999) [arXiv:hep- 
ph/9903358]. 

[18] A. M. Stasto, K. J. Golec-Biernat and J. Kwiecinski, Phys. Rev. Lett. 86 
(2001) 596 [arXiv:hep-ph/0007192]. 

[19] I. Balitsky, Nucl. Phys. B 463, 99 (1996). 

[20] Yu.V. Kovchegov, Phys. Rev. D 61, 074018 (2000). 

[21] J. Jalilian-Marian, A. Kovner, A. Leonidov, H. Weigcrt, Nucl. Phys. B 504, 
415 (1997); ibid. Phys. Rev. D 59, 014014 (1999); ibid. 034007 (1999), 
Erratum, ibid. 099903 (1999); J. Jalihan-Marian, A. Kovner, H. Weigert, Phys. 
Rev. D 59, 014015 (1999); A. Kovner, G. Milhano, H. Weigert, Phys. Rev. D 
62, 114005 (2000); H. Weigert, Nucl. Phys. A 703, 823 (2002); E. lancu, A. 
Leonidov, L.D. McLerran, Nucl. Phys. A 692, 583 (2001); ibid. Phys. Lett. 
B 510, 133 (2001); E. Ferreiro, E. lancu, A. Leonidov, L.D. McLerran, Nucl. 
Phys. A 703, 489 (2002). 

[22] A. H. Mueller, Nucl. Phys. B415 (1994) 373-385. 

[23] S. Munier, R. B. Peschanski, Phys. Rev. Lett. 91 (2003) 232001. [hep- 
ph/0309177]. 

[24] L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994); Phys. 
Rev. D 49, 3352 (1994). 

[25] J. L. Albacete, N. Armesto, J. G. Milhano et al, Phys. Rev. D71 (2005) 
014003. [hep-ph/0408216]. 

[26] J. L. Albacete, N. Armesto, A. Kovner et al, Phys. Rev. Lett. 92 (2004) 
082001. [hep-ph/0307179]. 

[27] A. H. Mueller, Nucl. Phys. B 643 (2002) 501 [arXiv:hep-ph/0206216]. 

[28] D. N. Triantafyllopoulos, Nucl. Phys. B648 (2003) 293-316. [hep-ph/0209121]. 

[29] J. P. Blaizot, F. Gelis and R. Venugopalan, Nucl. Phys. A 743 (2004) 13 
[arXiv:hep-ph/0402256]; ibid. 57 [arXiv:hep-ph/0402257]. 



26 



[30] J. Jalilian-Marian, Y. V. Kovchegov, Phys. Rev. D70 (2004) 114017. [hep- 
ph/0405266]. 

[31] C. Marquet, Nucl. Phys. A 796 (2007) 41 [arXiv:0708.0231 [hep-ph]]. 

[32] E. lancu, K. Itakura, L. McLerran, Nucl. Phys. A708 (2002) 327-352. [hep- 
ph/0203137]; A. H. Mueller, D. N. Triantafyllopoulos, Nucl. Phys. B640 
(2002) 331-350. [hep-ph/0205167]. 

[33] F. Caola and S. Forte, Phys. Rev. Lett. 101 (2008) 022001 [arXiv:0802.1878 
[hep-ph]]. 

[34] D. Kharzeev, E. Levin and M. Nardi, Phys. Rev. C 71 (2005) 054903 
[arXiv:hep-ph/0111315]. 

[35] J. L. Albacete, Phys. Rev. Lett. 99 (2007) 262301 [arXiv:0707.2545 [hep-ph]]. 

[36] J. P. Blaizot, T. Lappi and Y. Mehtar-Tani, Nucl. Phys. A 846, 63 (2010) 
[arXiv:1005.0955 [hep-ph]]. 



27 



